Acoustic attenuation probe for fermion superfluidity in ultracold atom gases 
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Dilute gas Bose- Einstein condensates (BEC's), currently used to cool fermionic atoms in atom 
traps, can also probe the superfluidity of these fermions. The damping rate of BEC-acoustic excita- 
tions (phonon modes), measured in the middle of the trap as a function of the phonon momentum, 
yields an unambiguous signature of BCS-like superfluidity, provides a measurement of the super- 
fluid gap parameter and gives an estimate of the size of the Cooper-pairs in the BEC-BCS crossover 
regime. We also predict kinks in the momentum dependence of the damping rate which can reveal 
detailed information about the fermion quasi-particle dispersion relation. 



5h 
43 



C 

o 

> 
m 

(N 

in 
o 
-i— > 



T3 

c 

o 
o 



X 
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Recently, cold atom physics entered [l| the previously 
inaccessible Bose-Einstein condensate (BEC) regime 
of fermion superfluidity. A magnetic field sweep through 
a Feshbach resonance induced the fermions to pair up into 
diatomic molecules which subsequently occupied a BEC 
(i.e. forming pairs hardly distinguishable from diatomic 
molecules - Shafroth pairs) . Current experiments aim at 
taking the fermion supcrfluid into the Bardeen-Cooper- 
Schrieffer (BCS) limit of ordinary superconductors, in 
which the fermion pairs are larger than the average inter- 
particle distance (Cooper pairs). The crossover region 
near the detuning of diverging scattering length, is of 
fundamental interest 0]. The cold atom crossover ex- 
periments [H| illustrate the need for unambiguous probes 
of fermion superfluidity. Sympathetic cooling immerses 
the fermionic atoms in colder bosonic atoms, often occu- 
pying a BEC. The system then resembles the condensed 
3 He- 4 He liquid mixtures In contrast to the helium 
liquid mixtures in which strong interaction effects drasti- 
cally reduce the boson-mediated fermion-fermion inter- 
actions Q, the cold atom mixture can exhibit BEC- 
mcdiated fermion superfluidity at attainable tempera- 
tures [|. 

In this paper, we show that the BEC simultaneously 
trapped with the fermion atoms, can provide a power- 
ful probe of fermion superfluidity. We study the damp- 
ing rate, r(g), of acoustic BEC-excitations (phonons) in 
the mixture as a function of the phonon momentum q. 
Acoustic attenuation 0] provides a signal for fermion 
superfluidity, determines the superfluid gap parameter, 
tests the validity of the weakly interacting quasi-particle 
description [l(J, estimates the size of the Cooper-pairs 
and probes the quasi-particle dispersion relation. The 
latter could be inferred from kinks which we predict to 
appear as an imprint of kinematical constraints. 
In linear response theory T(q) is related to density- 
density correlations. For a sufficiently large trap, the 
damping rate measured near the middle of the trap, 
where the particle densities are approximately constant, 
equals that of a BEC-phonon in a homogeneous mixture. 
The BEC-phonon is a Bogoliubov excitation with en- 



ergy Eg = hqcy/l + (£g) 2 , where c denotes the BEC 
sound velocity, c = ?t/(2ms^), £ is the BEC-healing 
length, £ = {XQupbcib) 1 ^ 2 , and as is the boson-boson 
scattering length. Using box normalization, the boson 
density ps corresponds to Nb bosons condensed in a 
volume fl, pb = Ng/fl. The system also includes a 
homogeneous mixture of spin + and — fermions, with 
p + = p_ = pp. Each single-component fermion system 
has Fermi wavenumber kp = (6n 2 pp) 1 ^ 3 , Fermi-velocity 
vp = hkp/mp and Fermi energy ep — h 2 k F /(2mp). 
The boson-fermion interaction is spin dependent, giv- 
ing different boson-fermion scattering lengths, a+ and 
a_, and a fermion-boson interaction Hamiltonian Hi„t = 
{2-Kh 2 a)(l/m B + l/mp)^ 1 £) k ,• P-k^jPfw wner e j = 



represents the boson (fermion) density opera- 
a/aj. 



tor, a = \J~o?- + arid aj = a/cij. We use perturba- 
tion theory around noninteracting fermion-boson prod- 
uct states, \)p (8 \)b- The _ffj„ t -interaction then anni- 
hilates a BEC-phonon |q) b and creates a fermion exci- 
tation that is a collective mode or a quasi-particle pair 
of opposite spins. By ignoring initial states with mul- 
tiple BEC excitations, we assume sufficiently low tem- 
peratures. By ignoring final states with two (or more) 
phonons, we make a relative error oc (<7£) 3 (as/£)(8/37r), 
but we take (<?£) < (C/a B ) 1/3 . Using B(0|p_ q |q) B = 



N B {h 2 q 2 /2m B )/(huj q 3 ), T(q) takes the form T(q) 
a*OjT\j(q) with 



r«(q) = 



8 (r /4^ g 

3 y/i + m 2 



PF 



dt 
h/e F 



(1) 



x exp[-i(qT-wfi)] (pf(r,t)pf (0,0)) 

where r = 47ra 2 p B w F [(l + mp/mB)/2] 2 , proportional 
to the rate at which a particle of velocity vp scat- 
ters off hard-core targets of radius a, distributed at 
the BEC density. By varying a + or a_, one can, in 

principle, distinguish the r++ and contributions. 

In a Cooper-paired fluid, the (p+pZ) correlations de- 
pend on the microscopic superfluid structure, embodied, 
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for instance, by the order parameter (i.e. the expec- 
tation value of annihilation field pairs, (t/i+V-) 7^ 0). 
The equal time ±-density fluctuation, (p F (r)p F (0)) ~ 
($ + (r)$(0))($l(r)$t(o)) = F(t)F*(-t), where F(r) = 

(0)) plays the role of fcrmion pair wavefunc- 
tion. Equation |T]) also suggests that one can probe di- 
rectional information if and when cold atom experiments 
realize anisotropic fermion superfluidity (such as p-wave, 
d-wave or LOFF-pairing 

HH El)- To make the quasi- 
particle pair explicit, we insert a summation over such 
states in the density-density correlation function. The 
temporal integral gives a delta function that expresses 
the kinematics of the pair creation process, 



r«(q) 



(8r /3)£g ep_ 

V 1 + (£<?) 2 p p 



d 3 k 



(2) 



where E F denotes the quasi-particle dispersion relation 
and gij (k; q) is the weight for the creation of quasi- 
particles of momenta k and q— k. The generalized Wick's 

theorem leads to g ++ = g__ = Hk(l — ?tq-k)j an d g-\ = 

g.+ = F^F* with F(r) = (2tt)- 3 / d 3 fc^ k exp(zk ■ r), 



and nk 



c j,k}- Here, c(c^) denote fcrmion annihila- 
tion (creation) operators. 
The form of Eq. © is quite general but the lack of corre- 
lations resulting from Wick's rule, is not. We can under- 
stand the role of F as pair wavefunction by comparing 
Eq. |(2J) to the Fermi Golden-rule rate at which a par- 
ticle of momentum q and energy scatters to zero 
momentum by breaking the bond of stationary diatomic 
molecules of (+) and (— ) atoms. In the initial molecules, 
these atoms occupy bound vibrational states with wave- 
function ifk and binding energy E\y. The interaction that 
breaks the molecular bond, Hi„ t , creates a pair of free 
atoms, the (+) atom with momentum k and energy efc, 
and the (— ) atom with momentum q — k and energy 
e |q-k|- The corresponding interaction matrix element 
is proportional to a-tp^ + a+ Pk-q, so that the size- 
dependent contribution to the rate is 

— cx a + a- [ d 3 k(p k (pl._ cl 5(E (l -E b -e k - e|k-q|) , (3) 

Tq J 



which has the form of Eq. @ if we associate F^ — > ip^ 

and E F — > Eb/2 + e^. Similar to the g_| contribution 

to Eq. ((2]), Eq. ([3]) stems from an interference of two pro- 
cesses in reaching the final state. In the first process, 
atom (+) with momentum k and atom (— ) with momen- 
tum — k in the molecule reach the final state because the 
(— ) atom receives a momentum kick q. In the second 
process, the (+) and (— ) atoms in the molecule have ini- 
tially momenta k — q and q — k, respectively, and the 
(+) atom receives a momentum kick q. 
The role of the fcrmion density of states, D(E) = 
(2ir 2 )~ 1 k 2 (E)/\dE F /dk\, becomes apparent when we 



convert Eq. ([2]) to an integral over quasi-particle energy. 
Targeting s-wave pairing, we assume a spherically sym- 
metric pairing wavefunction, _F q = F q . Introducing the 
angle 9 between k and q, we write 



E F (k) 



Vk 2 



2kqcos9 



|q — k\j 5(cos9 — cos 6j) 
k 3 q\dE(K)/d K \ K ^_^ 



(4) 



where the j-subscript enumerates the energy and mo- 
mentum conserving quasi-particle pairs. Whether or not 
the fermion dispersion relation supports a local minimum 
(i.e., whether or not the fermion supcrfluid resides in the 
BCS regime), affects the number of such pairs. If E F 
has a local minimum, we refer to it as the "gap" A, 
E F — A and we refer to the corresponding surface 
in momentum space, k = k m , as the gap surface. We 
can then assign two momentum magnitudes for the same 
quasi-particle energy E: k < (E) inside the gap surface 
(k < (E) < k m ) and k>(E) outside (k>(E) > k m ). The 
phonon annihilation can generate four quasi-particle pair 
types, j = (<, <), (<, >), (>, <), (>, >), where we indi- 
cate the momentum position of the (+) particle first, and 
that of the (-) particle second. Depending on the energy 
distribution, some pair creation processes j arc excluded 
kinematically. We can picture the constraints geomet- 
rically for the j-pair creation by considering the spher- 
ical surfaces in momentum space of radii kj t +(E) and 
kj-(E'), with E+E' = E^ . As the (+) momentum k lies 
on the first surface, the ( — ) momentum q — k on the sec- 
ond, the surfaces must be connected by the q-momcntum, 
requiring k j>+ (E) + %,_(£') <q< \k jt+ (E) - kj-{E% 
Then, the damping rates read 



7+~(g) 

\ 3PF J 

7++(g) 

( 87T^e F \ 
V 3PF J 
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D 3 {E)F {E) Dj(E')Fj(E') 



k 3 {E>) 



(5) 



D J (E)n 3 (E)D 3 (E')[l~n 3 (E')} 



kj(E) 



where ^(q) = [Ty («) y/1 + (£g) 2 ]/r . What is the 
information content of r(g)? The gap prevents soft 
phonons with E? < 2A from damping through pair 
excitation. If these soft phonons are kinematically 
prohibited from exciting collective fcrmion modes be- 
cause the latter's velocity (~ vp for the Anderson- 
Goldstone mode in the BCS limit) exceeds c, then 
they can only decay through Bcliaev damping which 
gives a negligible damping rate at long wavelengths 
(r Bc iiacv(g) ~ (6 v / 7r/5) v / ^a|'(g'c)(^) 4 ). More ener- 
getic BEC phonons damp by radiating quasi-particle 
pairs, and the higher the density of states, the more an 
energy interval contributes to the damping. At the gap 
surface the density of states is infinite, although the di- 
vergence is sufficiently slow (cx y E 2 — A 2 ) to ensure that 



3 



- 

a 
g 
IS 




1.5 



FIG. 1: Scaled damping rates components, 7++ and 7--, 
in the BCS limit, as a function of the phonon energy, Eg , 
measured in units of A, the superfluid gap parameter. Here, 
we have 7^ = Fyy/l + (£g) 2 /To- 



the integral |(SJ) remains finite. The divergence plays a 
prominent role when both quasi-particles form on or near 
the gap surface, corresponding to a phonon energy that 
slightly exceeds 2A. While the damping rate remains 
finite, the diverging density of states causes a discon- 
tinuous jump of the damping rate at the (edge) phonon 
momentum q& for which E^ & = 2A, see Fig. Q] Con- 
versely, the discontinuity in the damping rate points to- 
wards a diverging density of states, providing a signature 
of fermion superfluidity on the BCS side of the crossover. 
This sharp feature also provides a sensitive way of deter- 
mining the superfluid gap from the edge momentum q&, 
A = (AbPb)(£<7a)\/1 + (<7a£) 2 - O n the other hand, the 
observation of a continuous onset of the damping with in- 
creasing q does not necessarily imply that the fermions, 
if superfluid, have entered the BEC-side of the crossover. 
Indeed, as the system nears the crossover from the BCS 
side in an adiabatic sweep through a Feshbach resonance, 
we expect k m to continuously tend to zero, whereas cja 
monotonically increases. Beyond the detuning for which 
k m = ?A, 9 A is too large to fit inside the gap surface. In 
that case, the creation of both quasi-particles on the gap 
surface, the event that leads to the discontinuous jump in 
damping rate, is kincmatically prohibited. The depen- 
dence of the damping rate on Cooper-pair size, suggested 
by the analogy with breaking the stationary molecule 
bond, Eq. ((3]), becomes apparent when we expand the 
damping rates near the edge. We make two assumptions: 
(i) (E^) 2 « A 2 + h 2 v 2 (k — fc m ) 2 , in which v characterizes 
the curvature of the dispersion near the gap. At k = k m , 
F(k m ) = F m and n(k m ) = n m . (ii) The fermion system is 
sufficiently far from the crossover point to ensure that all 
four quasi-particlc pair events arc kincmatically allowed. 



Therefore, the damping rates near q 
-(?) 
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2^k F ){A/[m F v 2 /2])F 2 l 



<7A, (q > <?a) arc 
(6) 
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2n 2 (£k F )(A/[m F v 2 /2])n m (l 
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2A 



A 



- «m) 

1 — 2n r , 



(7) 



(Hv) 2 1 — n n 



"A- 



k=k„ 



Note that the ratio of the edge values of the (+, +) 
and (+, — ) damping rates connect n m and F m : 
r+-(<?A)/r ++ (qA) = n m (l - n m )/F 2 n . This ratio equals 
one in a mean-field description for short-range pairing 
interactions because then F m — n m = 1/2. This ra- 
tio can deviate significantly from one when long-range 
interactions are present, particularly in the crossover re- 
gion. If F m = n m = 1/2 and the (+,+) and (+, — ) 
edge damping rates are identical, then one can mea- 
sure v from their absolute value, since these are pro- 
portional to (mv 2 /2A) and A follows from q&- From 
the slope of the scaled damping rate 7_|_ > _(g), set out as 
a function of the phonon energy E^ , one can extract 
the square of the coherence length, — V 2 z F/ !; /F) <; \k=k m 
I d 3 rF(r)r 2 exp(ik m • r)/ J d 3 rF(r) cxp(ik m • r), where 
|k m | = k m . Near crossover, where the extent of F 
drops below k^ 1 , this coherence length tends to the 
range of F. In the BCS-limit, k m — > k F . v — > v F and 
— W 2 Fk / Fk\k=k m — * (hv F /A) 2 so that the slope tends to 
— 1/4, as shown in Fig. [TJ Near the crossover q& grows 
and, therefore, it is likely that the kincmatical condi- 
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FIG. 2: Near-edge pair breaking contribution to the scaled 
damping rate 7 (7 = r^/1 + (£g) 2 /ro), as a function of the 
phonon momentum, q. Here we assumed a+ = a_, c/vf = 
0.5 and varied the fermion-fermion scattering lengths bf or 
77 = (/cf^f) -1 - The edge momentum gives the superfluid gap, 
the slope near the edge relates to the size of the Cooper-pairs 
and the discontinuities in the derivatives (kinks) are imprints 
of the kinematical constraints of the pair breaking process. 
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tions \k jl+ (E) - kj-{E')\ < q < k jt+ (E) + kj-(E') are 
violated near the edge. At the minimal phonon energy 
Eq for which quasi-particle pairs j are kinematically ex- 
cluded, the damping curve exhibits a kink. We expect 
the kinks to be smoothed out if the quasi-particles have 
finite lifetimes. In Fig. [21 we show T(q) for a zero temper- 
ature Cooper-paired system at different detuning values 
during an adiabatic sweep through the resonance. The 
curves were calculated from a Keyldish-like wavefunction, 

_ J, t t 

e °L t <Pfc c k,+ c -k,- 1 0) , where we choose (f> k as a Yukawa 
wavefunction with the range given by the BCS coherence 
length, R = A/hvp. While the assumed wavefunction 
cannot describe the system quantitatively, the resulting 
damping rates illustrate the relevant physics. 
To demonstrate the practicality of cold atom acoustic at- 
tenuation, we consider a realistic two-component fermion 
( 6 Li)- BEC ( 23 Na) mixture (a B = 3.4nm,£ = 0.14/zm), 
at temperature T = 50 nK, trapped by harmonic oscil- 
lator potentials within a radius Rq = 50 /xm with peak 
densities pp = 10 12 cm -3 and pb — 3 x 10 14 cm~ 3 . We 
assume that a = 5nm and that the system has a gap 
A ~ ep/5. For these parameters, XbPb = 266nK, 
ep = 610 nK (A = 122 nK), the relevant velocities are 
vp = 4 cm/sec and c = 1 cm/sec, the time scale for BEC- 
phonon scattering is T^" 1 = 0.17 msec during which time 
the BEC-excitation travels a distance T^c ~ 1.7 pm, 
and the critical phonon momentum for pair breaking, q&, 
(obtained from (g A £) 2 = -1/2 + y/1/4 + [A/(X B p B )} 2 ~ 
0.178) is equal to q& — 3 pm^ 1 . Creating the wavepackct 
by overlapping two laserbeams over a region L of 5 pm 
gives a relative momentum spread Sq/q& ~ l/(Lq&) ~ 
0.067. Measuring the phonon damping or scattering af- 
ter a time less than or equal to 0.5 msec, this wavepacket 
would have traveled over less than 5 pro., so that the mea- 
surement can be carried out within a radius L mc&s = 
5 pm [l]|. The relative density variations, only one per- 
cent for L meas /Ro ~ 0.1 since 5p B /pB ~ (L mcas / 'R ) 2 , 
Spp/pF ~ (3/2)(L moas /i? ) 2 , smooth out the sharp fea- 
tures of the damping curve but only over a few precent 
of the phonon momentum, less than the relative momen- 
tum uncertainty due to the finite size of the wavepackct 
~ 5 — 7%. Even with a smoothing of ten percent of the 
momentum, we expect the most abrupt kinks shown in 
Fig. [2] to be clearly visible. At finite temperature, ther- 
mal scattering contributes a smooth background that is 
considerably smaller than the pair breaking term at suf- 
ficiently low temperatures, such as in the above example. 
At A < k B T, the thermal occupation of fermion quasi- 
particles is exponentially suppressed and so is the cor- 
responding contribution to the phonon scattering. The 
scattering by thermally excited BEC phonons is not sup- 
pressed, but still small: estimating its magnitude roughly 
as 7thcrm ~ 4ira 2 p B c B f where / is the fraction of ther- 
mally depleted bosons, / - 0.02 (from / - (T/T BE c) 3/2 , 



where the BEC critical temperature, Tbec ~ 20 pK at 
p B = 3 x 10 14 cm" 3 ) gives 7thcrm ~ 0.03r . 
In summary, we have shown that the damping rate of an 
acoustic excitation of an overlapping BEC can probe the 
superfluidity of an ultra-cold fermionic gas. The momen- 
tum dependence of this rate yields an unambiguous sig- 
nature of BCS-like superfluidity, allows a measurement 
of the superfluid gap parameter and of the correlation 
length in the BEC-BCS crossover regime. Kinks in the 
momentum dependence of the damping rate will reveal 
information about the quasi-particle dispersion relations. 
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